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We present 2D simulations of the cooling of neutron stars with strong magnetic fields 
Oh' (B > 10 13 G). We solve the diffusion equation in axial symmetry including the state 

of the art microphysics that controls the cooling such as slow/fast neutrino processes, 
superfluidity, as well as possible heating mechanisms. We study how the cooling curves 
depend on the the magnetic field strength and geometry. Special attention is given to 
discuss the influence of magnetic field decay. We show that Joule heating effects are very 
large and in some cases control the thermal evolution. We characterize the temperature 
anisotropy induced by the magnetic field for the early and late stages of the evolution of 
isolated neutron stars. 
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1. Introduction 

The observed thermal emission of neutron stars (NSs) can provide information 
about the matter in their interior. Comparing the theoretical cooling curves with 
observational data 1,2 one can infer not only the physical conditions of the outer 
region (atmosphere) where the spectrum is formed but also of the poorly known 
interior (crust, core) where high densities are expected. 

There is increasing evidence that most of nearby NSs whose thermal emission is 



visible in the X-ray band have a non uniform temperature distribution 3 ' 4 . There is 
a mismatch between the extrapolation to low energy of the fits to X-ray spectra, and 
the observed Rayleigh Jeans tail in the optical band (optical excess flux), that cannot 
be addressed with a unique temperature (e.g. RX J1856.5— 3754 5 , RBS 1223 6 , 
and RX J0720.4-3125 7 ). 

A non uniform temperature distribution may be produced not only in the low 
density regions 8 , but also in intermediate density regions, such as the solid crust. 
Recently, it has been proposed that crustal confined magnetic fields with strengths 
larger than 10 13 G could be responsible for the surface thermal anisotropy 9,10 . In 
the crust, the magnetic field limits the movement of electrons (main responsible 
for the heat transport) in the direction perpendicular to the field and the thermal 
conductivity in this direction is highly suppressed, while remains almost unaffected 
along the field lines. 
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Moreover, the observational fact that most thermally emitting isolated NSs have 
magnetic fields larger than 10 13 G implies that a realistic cooling model must include 
magnetic field effects. In a recent work 11 , first 2D simulations of the cooling of 
magnetized NSs have been presented. In particular, it has been stated that magnetic 
field decay, as a heat source, could strongly affect the thermal evolution and the 
observations should be reinterpreted in the light of these new results. We present 
the main conclusions of this work next. 

2. Non-uniform temperature distribution induced by magnetic 
field 

We consider two baseline models 11 : Model A, a low mass NS with M — 1.35 M Q and 
Model B, a high mass NS with M — 1.63 M Q . These two models correspond to the 
minimal cooling scenario (controlled by modified Urea neutrino emission) and the 
fast cooling scenario (where direct Urea operates), respectively. We consider crustal 
confined magnetic fields, keeping the geometry fixed and varying the magnetic field 
strength at the pole (B). 




Fig. 1. Cooling of strongly magnetized NSs for the Model A. Temperatures TJ, (left figure) and T s 
(right figure) at the pole and at the equator as a function of the age t (left panels). Corresponding 
evolution of T /T v ° le vs 9 in the right panels. 



For a given magnetic field, we solve the diffusion equation in axial symmetry, con- 
sidering an anisotropic thermal conductivity tensor k. The ratio of its components 
along and perpendicular to the field can be defined in terms of the magnetization 
parameter (lubt) as kI/k^t = 1 + (wbt) 2 , where r is the electron relaxation time 
and wb is the electron cyclotron frequency. When u>bt 3> 1 the magnetic field ef- 
fects on the transport properties are crucial. We show the results in Fig. Q] (on the 
left), where the temperature at the base of the envelope (at 10 9 g/cm 3 ) is shown 
as a function of the age t. We see an inverted temperature distribution with cooler 
polar caps and a warmer equatorial belt at t < 500 yr for the Model A. Similar 
qualitative results are found for Model B. 
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In Fig. [T] (on the right) we show for Model A the corresponding cooling curves 
but for the surface temperature T s . We see that the anisotropy found at the level of 
Tf, does not automatically result in a similar T s distribution: the blanketing effect of 
the envelope overrides the inverted temperature distribution found at intermediate 
ages. Thus, the equator remains always cooler than the pole, and only at early times 
and for strong fields we find larger surface temperatures in middle latitude regions. 

3. Magnetic field decay and Joule heating 

For the magnetic field decay, we assume the approximate solution of the diffusion 
equation 

p = p CX P M/ T Ohm) m 

o 1 + ^ (1 _ exp( _ Vrohm)) y ) 

where 7ohm is the Ohmic characteristic time, and the typical timescale of the fast, 
initial stage is defined by Tall- 
in the cooling curves, there is a huge effect due to the decay of such a large 
field: as a consequence of the heat released, T s remains much higher than in the 
case of constant field (Fig. ^ . The strong effect of the field decay is evident for all 
of the pairs of parameters chosen. Notice that T s of the initial plateau is higher for 
shorter THaii , but the duration of this stage with nearly constant temperature is also 
shorter. After t ~ t k&ii , there is a drop in T s due to the transition from the fast 
Hall stage to the slower Ohmic decay. 




Fig. 2. Cooling of strongly magnetized NSs with Joule heating with Bo = 5 X 10 G for Model 
A (left panel) and Model B (right panel). We have set TQhm = 10 6 yr. 
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4. Conclusion: Towards a coupled magneto-thermal evolution 

The main result of this work is that, in magnetized NSs with B > 10 13 G, the decay 
of the magnetic field affects strongly their cooling. In particular, there is a huge 
effect of Joule heating on the thermal evolution. In NSs born as magnetars, this 
effect plays a key role in maintaining them warm for a long time. Moreover, it can 
also be important in high magnetic field radio pulsars and in radio-quiet isolated 
NSs. As a conclusion, the thermal and magnetic field evolution of a NS is at least 
a two parameter space (Fig [3]), and a first step towards a coupled magneto-thermal 
evolution has been given in this work. 



Magnetars | 




Fig. 3. Coupled magneto-thermal evolution of isolated neutron stars 12 : T s for the hot compo- 
nent as a function of B and t. Observations: squares for magnetars (B > 10 14 G), triangles for 
intermediate-field isolated NSs (10 13 G< B < 10 14 G) and circles for radio pulsars (B < 10 12 G). 
Corresponding cooling curves in solid, dashed and dotted lines, respectively. 
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